Asymmetry of the Hamiltonian and the Tolman's length 
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Using the canonical transformation of the order parameter which restores the Ising symmetry of 
the Hamiltonian we derive the expression for the Tolman length as a sum of two terms. One of 
them is the term generated by the fluctuations of the order parameter the other one is due to the 
" , entropy. The leading singular behavior of the Tolman length near the critical point is analyzed. The 

i-^. ■ obtained results are in correspondence with that of M.A. Anisimov, Phys. Rev. Lett., 98 035702 

(2007). 
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Introduction 

oo : 

As is known (see e.g. JU) the Tolman length can be defined as the correction caused by the fact that the equimolar 
f— - n , surface does not coincide with the surface of tension for a small droplet. The origin of such a difference is the 
' asymmetry of the phase coexistence in terms of the density variable [2] . In its turn it is due to the asymmetry of 
the Hamiltonian. In [2] it was shown that within the square gradient model the nonzcroth value of the Tolman 
length is caused by the asymmetry of the density functional (Helmholtz free energy). For the symmetrical models 
such as lattice gas the Tolman length vanishes identically. Such a situation resembles the one with the issue of the 
rectilinear diameter and its singularity . This effect is due to the asymmetry of the Hamiltonian as a functional 
of the order parameter such as density. The effects of asymmetry are consistently treated within the approach based 
on the canonical form of the Hamiltonian [J, H Q . In such an approach both the linear [H, 0] and nonlinear Q 
mixing of the thermodynamic fields are treated uniformly on the basis of the isomorphism principle. In particular, the 
effects caused by the asymmetry of the Hamiltonian are interpreted as the sequence of improper choice of the order 
parameter. For proper, canonical order parameter rj the Ising symmetry with respect to transformation rj — > — r; is 
restored. The main difference between "complete scaling" approach and the proposed approach the canonical form 
of the Hamiltonian lies in the fact that one does not need to use the three scaling fields but work directly with the 
Hamiltonian. In particular, it allows to make some conclusions about the amplitudes of the singularities. 

Recently, the question about critical behavior of Tolman length is analyzed in [l(| within the "complete scaling" 
approach. The leading singularities were obtained. 

In this paper we demonstrate that the canonical formalism leads to the essentially the same results obtained in 
in addition giving the possibility to relate the isomorphic variables with the microscopic Hamiltonian. Basing on the 
expressions for the Tolman length given in 0, [TTj] we obtain the expression for the Tolman length which explicitly 
C*~) . demonstrates the role of the asymmetry of the Hamiltonian. 

The structure of paper as follows. In Section Q] we give short outline of the procedure of the reduction of the 
Hamiltonian to the canonical form. In Section |TT] we use the canonical transformation for structuring the Tolman 
"xl" ' length. 
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I. THE REDUCTION OF THE HAMILTONIAN TO THE CANONICAL FORM 



We consider the case of 2-nd order phase transition. The typical examples are Ising model and simple molecular 
j_j liquids. For these systems the Hamiltonian takes a form 

H[p(r)} = H ql [p(r)] + J H loc {p{v))dv , (1) 

where 

oo 

#UMr);K}) = £% n (r). (2) 

ft 

n—l 

For convenience we include (3 — -^Kp into the Hamiltonian. In a case of simple liquids, for example, the coefficients a n 
are definite functions of the chemical potential fi (in fact the difference of the chemical potential and its value /i(T) 
at the coexistence curve) and the temperature T if the nontrivial reference system is used [13] • 
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Due to locality, for every point we can write: 

p(r) = F( V (r j) , (3) 
then for the integrand in the partition function of the system we have 

cxp (-<r%)) = / HP- Ffo))exp (-Ht n \p)) dp, (4) 

where 

#Zo C W = A iV + -y»7 + —V ■ (5) 

Here we give the procedure of transforming the Hamiltonian to canonical form using the ideas of the Catastrophe 
Theory [13 ]. The link of the CF with the CT is the transformation p — > rj(p) of the initial order parameter. Similar 
idea of transformation of the variable reducing the distribution to simpler (gaussian) form was used in [l4j . 

It is expedient to represent the local part of the fluctuational Hamiltonian as the sum of even "(+) " and odd "(— ) " 
parts: 



H loc (p(r)) = £%"(r) = H£J(p(r)) + h£{ P {t)) (6) 

n—1 

For the local part the Hamiltonian we can write: 

m = G[ P ] , (?) 

where 



F[fj)=j exp(-H% n \z))dz 


p 

G[p]= J cxp(-H loc (z))dz . 



The coefficients A\ and A2 of the canonical form are determined basing on the equalities 

F[+oo; Ai,A 2 , A A ] = G[+oo; p, T] 

F[-oo; A X ,A 2 , M] = G[-oo; p, T] (8) 

which provide the bijectivity of the transformation given by Eq. ([7])- To be specific we assumed that the coefficients 
of the initial Hamiltonian depend on the "laboratory" variables chemical potential p and the temperature T. Two 
conditions ijHJ) arc not sufficient to fix three coefficients Ai,A2,Ai as functions of the "laboratory" variables. The 
coefficient A4 can be fixed according to some additional condition since it is assumed that A4 7^ 0. 

Let us describe how the coefficients of the canonical form relate with the laboratory variables. Note that the 
Hamiltonian ^ is usually based on the mean-field equation of state. Therefore neglecting the fluctuations the 
coefficients of the canonical form can be found from the condition of invariance of the CP locus within the local 
(mean- field) approximation (the coefficient 03 also vanishes at the CP due to stability reason [IB]): 

A 1 {p c ,T c ) = 0, A 2 {p c ,T c ) = ai {p c ,T c )=0, a 2 {p c ,T c ) = 0. (9) 

These equations fix the value of A4: 



+00 



A 



4 4 



+00 

drf = J exp (-Hioc(p; 01 = 0, a 2 = 0)) dp , (10) 



which gives 



= r~4 • (ID 

+oo \ 

J exp(-H loc (p;ai = 0, a 2 = 0))dp\ 
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But bearing in mind the account of fluctuations we do not need such an approximation. Moreover the fluctuations 
change the locus of the critical point. Taking into account the subsequent renormalization of the Hamiltonian ([5]), it 
is natural to put A4 equal to the renormalized interaction constant u* of the nongaussian fixed point: 

At = u* (12) 

and use the representation A 2 = r* + r, where r* is the coordinate of the nongaussian fixed point r* = — | A 2 + o(e) 
and A is the momentum cutoff [16j . In particular the critical point is determined by: 

A x {n c ,T c ) =0, A 2 (fi c ,T c ) =r*, 

Since the transformation (U|) is smooth, the coefficients A^ are the smooth functions of the laboratory variables. In 
particular at the coexistence curve for liquid-vapor transition we have: 

A 2 (^T)\ Ai=0 = ar + o(r), (13) 

T -Tr 



a > 0, 



Thus the coefficients A\ and A 2 of the canonical form can be treated as the scaling fields. This is the approximation, 
because the higher gradient terms are omitted. But since such terms do not contribute to the renormalization of 
the local terms we can expect that such an approximation incorporates the main thermodynamic features of the 
interparticle interactions. Thus in variable rj the Hamiltonian functional takes the Landau-Ginzburg form: 

Hlg fo(r)] = J (A lV + ^ ry 2 + ^ ^ + \ (V ?? ) 2 ) dV , (14) 

where we scaled the square gradient term appropriately. Expanding Eq. (J7J) in a series we get: 

i/(r) = p(r) + \v 2 p(rf + ^T 3 p(rf + ^r 4 p(vf + . . . , (15) 

where all coefficients are functions of the "laboratory" variables (e.g. /x and T): 
T 2 =ai + Ax , 

r 3 =i(a 2 + ai 2 + A 2 +A\ + 1 ia 1 A 1 ) , (16) 

1117 3 7 

T4 = x a 3 + 77 a x a 2 + - + - A X A 2 + A^ + - a 2 A x + - a x 2 A x + 2 oi^i 2 + a x A 2 . 
3 2 6 6 2 6 

The approach proposed allows to treat correctly the effects caused by the asymmetry of the Hamiltonian. Note 
that all information about the asymmetry caused by the odd part of the local hamiltonian ti\ Q } is represented by the 
linear term of the local part of of the canonical form for the Hamiltonian with A\ 7^ and is encoded also into A 2 
and A4. 

The influence of asymmetry of the Hamiltonian is naturally to analyse via the representation: 

H loc {p{v)) = H^(p(r)) + H£(p(r)) , (17) 



^W) ^^-^ , (is) 

H (-) (p(r)) = HloMr)) +H l0c (-P(r)) (1Q) 

where superscripts (+) and (— ) stand for even and odd components of the function correspondingly. Assuming that 
the odd part of the Hamiltonian is "small" in comparison with the even part, from Eq. ||5J) we obtain: 

+00 

A ± = 1 J H£ (x) exp ( -H<£ {x))dx + o[ ) , (20) 
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where 



ci = 



■ erfc 



-4, 



(21) 



Neglecting the fluctuational shift of the critical point this result shows that the main contribution to A% comes from 
05 (assuming that the coefficients a2 n +i decrease with n). In [6| it was shown that the singularity of the rectilinear 
diameter is shared by both r 2 ' 3 and r 1_Q anomalies, which are generated by the asymmetrical part of the Hamiltonian. 
The result of [6J together with Eq. ([2D)) allows to relate the amplitudes of the r 2 ^ and r 1_a singularities with the 
specificity of the interparticlc interaction of the system. In particular in [6[ it was shown that these amplitudes 
have opposite signs. This is in correspondence with the results obtained in 0, From Eq. (|2"0]) it follows that 
perturbatively this sign is determined by the coefficient a§. The value T 2 can be considered as the asymmetry factor. 
Form Eqs. (fT6|) . (f2Tj) it follows that such a factor is determined by the asymmetrical part of the Hamiltonian. 
Within the same approximation the coefficient A 2 is determined implicitly as: 



2A A 



Ki 



8A 4 



exp (-#<+>(*) Jdx 



(22) 



where K q is the q-th order modified Bessel function of the 2-nd kind. 

The variables Ai and A 2 unify the description of different systems near the critical point. The relations Eqs. (I20p . 
(|22| relate the coefficients of the canonical form (|T4| with the laboratory variables. This dependence determines the 
critical amplitudes for the specific system. 

From Eq. (fT5")) it follows that for the average value of the density we can write: 



P V&q ' Vasym T" 



(23) 



where 



Veq = ( V (r))=±\A 2 fg v ^ 



A x 



\A 2 \P+"t 

is the equilibrium value of the canonical order parameter, and 



Wegner corrections 



IJasym 



-r 2 (r4 +s v ) , ^ = (^( r ))-( r7 ( r )) 2 = |A 2 | 1 -^^^ TT ) 



regular terms . 



(24) 



is the part of the order parameter generated by the asymmetry of the Hamiltonian. The representation (|23[) together 
with Eq. (|24|) gives the general basis for treating the asymptotical properties of the physical quantities caused by the 
asymmetry of the coexisting phases (e.g. liquid and vapor). In particular, Eq. (|24j) describes the rectilinear diameter 
singularity 

Also it leads to the corresponding representation of the dimensionless (isothermal) compressibility of the coexisting 
phases into symmetric and asymmetric parts: 



P 2 XT 



dp 
dp 



9Ai(p, T) 



dp 



d dA 2 (fi,T) 



8A 1 



dp 







dA 2 



( Veq Tjasym ) ~l~ * • • Xsym H~ Xasym 



(25) 



analogous to that obtained in 0] (see also [T(|). Despite similar ideology (the nonlinear transformation of the order 
parameter) the formalism of the canonical form of the Hamiltonian differs from the "complete scaling" approach 
within the "complete scaling" approach originally proposed in 0] to treat the singularity of the rectilinear diameter 
and to resolve the nature of Yang- Yang anomaly [18| ■ The proposed approach is based directly on the Hamiltonian. 
This leads to the prediction that both r^ q and s n contributions are generated by the asymmetry of the Hamiltonian 
and proportional to asymmetry factor T 2 - While in "complete scaling" approach they are in fact independent because 
of the phenomenological nature of the hypothesis of complete scaling [H, [l3|- ln addition the proposed approach 
predicts that these two contributions have opposite signs. This seems in correspondence with the estimates found 
[a lS( by processing both the experimental data and model systems. 

From Eq. (|25[) we see that in addition to the standard |t|~ 7 singularity which is the same for both coexisting phases 
we get the leading correction terms |r|^~ 7 and |r| 1_a_/3_7 which take opposite signs in these phases at the coexistence 
curve due to r] aS y m . This result coincides with that of [l(|. Since such terms according to [T3| determine the critical 
behavior of of the Tolman length we will consider the application of the canonical formalism and the representation 
(|2"3")l to this problem. 
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II. CRITIC ALITY OF THE TOLMAN LENGTH 



According to [2| the Tolman length as the coefficient of the asymptotic correction to the surface tension of a drop 
of radius R at R — » oo is as following: 

+00 +00 
J zp'(z)dz J zp' 2 {z)dz 

S ? = ^ ^ (28) 

/ P '(z)dz J fP{z)dz 
—00 — 00 

where p(z) is the equilibrium density profile of interphase coexistence. This profile is obtained basing on the mini- 
mization square gradient functional (see e.g. [l[): 

F[ P (r)}= J (y (Vp) 2 + /[p(r)]) dr. (27) 

which in fact is the LGH. In accordance with the result of previous section we perform the local canonical transfor- 
mation 

p(r)= V (r)-^T 2 r, 2 (v)+... , (28) 

where p = p(r)/p c — 1, which restore the symmetry of the functional (|27p in canonical variable r\. 

Note that according to the definition, the spatial profile of the canonical order parameter rj eq (z) for two phase 
coexistence is an odd function with respect to the interphase boundary which is defined as the "equi-77" surface: 

Veq(-Z) = ~Veq{z), (29) 

, ust like for any model with the even Landau- Ginsburg functional [2]. The approach based on the correlation functions 
19| gives the same result. The second term in Eq. (|23|) represents the asymmetry effects [HQ. In particular 
Eqs. (|2l?l) . J24p lead to the "r 2/3 -" and " T 1_Q -" anomalies of the the rectilinear diameter [f|. From Eq. (|2"3")l it follows 
that the phase coexistence profile of the density can be written as follows: 

P{z) = rjeq(z) + r] asym (z) + ... 

with obvious Ising like properties 

Veq(z) = -Veq(-Z) , Vasym(z) = Vasym{-z) , (30) 

because of the symmetrical form of the Hamiltonian in canonical variable rj. Substituting this expression into Eq. J26|) 
we obtain: 

Sp = ~r 2 Scan + ■ ■ • , (31) 

where 

Scan = S„ + S s . (32) 

Thus the amplitude value of the Tolman length is governed by the value of I^, which can be either positive or negative 
depending on the details of the microscopic interaction. Since ij eq oc |r|^ and rj asym oc \t\ 2 @ to the leading order we 
can write: 

I zd Vcq(z) J z r] eq (z)r]' e 2 q (z) dz J zds v [z) J z s' v (z) dr] eq (z) 

<L = -=^ 2^— , S s = -^ =- (33) 

/ drjeq{z) J V'e q (z)dz J dr] eq {z) J Ti' e 2 q (z)dz 

— oo — oo — oo — oo 

where r](z) and s v (z) are the equilibrium profiles of the canonical order parameter and the entropy correspondingly. 

Further we assume that the density profile varies over the correlation length £ as the only relevant characteristic 
spatial scale near the critical point. Below we give the ground to such an assumption using the rigorous thermodynamic 
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expression for 6 obtained in [Tl| ■ Then simple scaling consideration shows that the obtained contributions to Tolman 
length have the following leading singular behavior 



S v oc r^ v . 



5 S OC T 1 " a - P ~ l/ . 



This is exactly the result obtained in [T(|. The expression (|3"Tj) for the Tolman length allows to give the ground for 
the approximate expression 



Pd 



Ap 



(34) 



proposed in [T(| basing on the asymptotic critical behavior. Here pd — p '"' 2 ^ Pa '" is the rectilinear diameter and £ is 
the correlation length. From the point of view of Eq. (|23|) it is definitely right qualitatively. In order to obtain the 
critical amplitudes for the behavior of the 6 we use the rigorous thermodynamic expression for the Tolman length 
given earlier in 0, H3] and recently represented in "compressibility form " in : 



(BJL \ 
\ 3/i T J 



(Apr 

Here is the surface tension of planar interface. Substituting Eq 

d rj. 



(35) 



and Eq. into Eq. ([33)1 we obtain: 



Xas 



Pc {ApY 



So that to the leading order we have: 



s _ ^oo dAi 



dp, 



g„(0) 
5,(0) 



Pc 



d, 



4^2 



eq 



-1-0 



(36) 



(37) 



Both expressions for Tolman length Eq. (f26|) and Eq. (|35|) give the same asymptotic behavior provided that a oc 
oc |^4 2 | 2,y thus proving the assumption that the correlation length is the characteristic scale for the spatial profile 
of the density and the surface tension as the specific thermodynamic potential of the interphase surface. From Eq. (|37|) 
it follows that the amplitude of the Tolman length is determined by the asymmetry factor T2. 



Conclusion 

In present paper within the canonical approach (see 0,01) we show that the non zeroth value of the Tolman length is 
the effect of the asymmetry of the Hamiltonian of the system in density variable. Performing the transformation to the 
canonical order parameter for which symmetry of the Hamiltonian is restored we derive the invariant representation 
of the Tolman length in terms of the profiles of canonical order parameter r\[z) and canonical entropy s[z). Such 
a representation allows to analyse the asymptotic behavior of the Tolman length. The leading singular terms are 
generated by the two above mentioned contributions and proportional to oc t^^ v and oc T 1 - a -f i - v correspondingly. 
This is in correspondence with the results of 10] obtained within the "complete scaling" approach of Q. In fact the 
qualitative representation (|34|1 shows that the nature of the singularity of the Tolman length is determined by the 
singularity of the rectilinear diameter pd- 
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